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This  paper  presents  a  numerical  model  of  an  exhaust  heat  recovery  system  for  a  high 
temperature  polymer  electrolyte  membrane  fuel  cell  (HTPEMFC)  stack.  The  system  is 
designed  as  thermoelectric  generators  (TEGs)  sandwiched  in  the  walls  of  a  compact  plate- 
fin  heat  exchanger.  Its  model  is  based  on  a  finite-element  approach.  On  each  discretized 
segment,  fluid  properties,  heat  transfer  process  and  TEG  performance  are  locally  calculated 
for  higher  model  precision.  To  benefit  both  the  system  design  and  fabrication,  the  way  to 
model  TEG  modules  is  herein  reconsidered;  a  database  of  commercialized  compact  plate- 
fin  heat  exchangers  is  adopted.  Then  the  model  is  validated  against  experimental  data  and 
the  main  variables  are  identified  by  means  of  a  sensitivity  analysis.  Finally,  the  system 
configuration  is  optimized  for  recovering  heat  from  the  exhaust  gas.  The  results  exhibit  the 
crucial  importance  of  the  model  accuracy  and  the  optimization  on  system  configuration. 
Future  studies  will  concentrate  on  heat  exchanger  structures. 

Copyright  ©  2012,  Hydrogen  Energy  Publications,  LLC.  Published  by  Elsevier  Ltd.  All  rights 

reserved. 


1.  Introduction 

With  the  diminishing  reserve  of  fossil  fuels  and  the  rising 
concern  on  global  green-house  gas  emissions,  more  attention 
is  being  drawn  by  fuel  cells  [1].  Among  other  fuel  cell  tech¬ 
nologies,  high  temperature  polymer  electrolyte  membrane 
fuel  cell  (HTPEMFC)  power  systems  are  considered  to  be 
promising  in  practice  [2].  Due  to  the  high  operating  tempera¬ 
tures  of  HTPEMFCs  (160°C-180  °C),  the  waste  heat  is  of  high 
quality  and  thermoelectric  generators  (TEGs)  are  introduced 
herein  to  recover  this  exhaust  heat  for  electricity.  As  electric 
power  generators,  TEGs  are  advantageous  in  low  grade  heat 
recovery  [3].  As  opposed  to  traditional  energy  conversion 
devices,  they  have  very  simple  structures  and  no  moving 
parts.  In  addition,  they  are  environmentally  friendly,  silent, 
reliable,  exhibit  excellent  scalability  and  pure  DC  power 


output.  Currently,  TEGs  are  restricted  to  low  heat-to-power 
conversion  ratios.  They  are  primarily  competitive  in  niche 
applications,  such  as  exhaust  heat  recovery,  remote  power 
supply  and  back-up  power  units  [4,5].  If  their  efficiency  can  be 
improved  in  the  future,  for  example  from  the  present  device 
figure  of  merit  (ZT)  ~  1  to  ZT  >  2,  many  applications  will 
materialize  [6,7].  More  widespread  use  of  thermoelectrics 
requires  not  only  novel  materials  with  improved  performance 
characteristics  [8—10],  but  also  TEG  systems  properly  con¬ 
structed  [11-13]. 

Research  on  TEG  systems  can  be  categorized  into  optimi¬ 
zation  of  the  TEG  modules  themselves,  and  assembly  of  TEG 
modules  with  heat  transfer  structures.  There  are  numerous 
publications  considering  the  inner  structure  of  TEG  modules 
[14—17].  However,  as  mentioned  by  Hendricks  and  Lustbader 
[18],  the  whole  TEG  system  should  be  studied  as  an  integrated 
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solution,  because  TE  material  properties  and  heat  exchanging 
performance  are  tightly  intertwined.  Modeling  and  simulation 
plays  the  main  role  in  the  studies  referenced  in  the  next  section. 

Much  effort  has  been  made  to  study  a  single-/multi- 
element  TEG  module  with  its  external  heat  reservoirs  [19—21]. 
However,  in  practice,  a  concrete  TEG  system  is  always 
a  combination  of  TEG  modules  with  proper  heat  exchanging 
structures  in  very  large  scale  installations  to  duly  account  for 
the  magnitude  of  the  exhaust  heat  sources.  Suzuki  and 
Tanaka  compared  the  performance  of  15  flat  thermoelectric 
panels  [22]  and  6  different  geometrical  arrangements  of 
cylindrical  thermoelectric  tubes  [23].  Crane  and  Jackson  [24] 
proposed  a  numerical  model  of  a  practical  TEG  waste  heat 
recovery  system  with  an  advanced  cross-flow  heat  exchanger 
and  later  the  model  was  validated  in  experiments.  Based  on 
the  aforementioned  research  work,  Yu  and  Zhao  [25],  Smith 
[26]  and  Espinosa  et  al.  [27]  used  different  numerical  models  in 
a  TEG  system  investigation  for  different  waste  heat  sources. 
Recovering  exhaust  heat  from  fuel  cells  is  a  relatively  novel 
application  of  TEGs.  Chen  et  al.  [28]  proposed  a  module-level 
three  dimensional  TEG  model  in  ANSYS  FLUENT  for  the 
convenience  of  the  co-design  and  co-optimization  of  a  PEMFC- 
TEG  system.  The  PEMFC-TEG  system  was  also  preliminarily 
demonstrated  in  the  lab.  A  solid  oxide  fuel  cell-TEG  (SOFC- 
TEG)  hybrid  system  was  optimized  on  the  main  variables  and 
the  operating  conditions  in  [29]  using  a  zero-dimensional 
system  model.  Overall,  the  previous  work  has  significant 
attributes  in  analyzing  TEGs  in  practical  applications.  Some 
features  can  be  further  enhanced,  including  the  following: 

1)  Fluid  properties,  flow  condition,  heat  transfer  [20]  and  TEG 
performance  are  locally  resolved  on  each  segment  in  a  dis¬ 
cretized  TEG  heat  recovery  system  model  for  higher  model 
accuracy. 

2)  The  recovery  system  model  is  engineering-oriented. 
Factors  affecting  the  system  performance  are  able  to  be 
co-optimized  in  the  model. 

3)  Mature  heat  exchanger  structures  and  commercially 
available  TEG  modules  are  used  in  the  model  to  benefit  the 
system  optimization  and  fabrication. 

In  the  present  paper,  a  TEG  exhaust  heat  recovery  system  is 
modeled  and  optimized  for  the  cathode  exhaust  gas  of 
a  HTPEMFC  stack.  Basically,  it  is  a  system  similar  to  compact 
plate-fin  heat  exchangers  but  with  TEG  modules  mounted  on  the 
walls.  Cathode  exhaust  gas  and  coolant  water  are  the  working 
media.  A  numerical  model  for  the  system  is  proposed  and  vali¬ 
dated.  The  model  is  intended  to  optimize  the  system  design  and 
operation  with  higher  accuracy  and  flexibility  as  well  as  to  assist 
the  system  fabrication.  Sensitivity  analyses  are  then  carried  out 
to  identify  the  parameter  priority  for  later  configuration  opti¬ 
mization.  Finally,  the  system  configuration  is  optimized. 


2.  Model  description 

The  chosen  architecture  of  the  TEG  exhaust  heat  recovery 
system  is  illustrated  in  Fig.  1.  The  compact  heat  exchanger 
housing  is  placed  in  the  middle  for  the  exhaust  gas.  Two 
aluminum  blocks  are  settled  at  the  top  and  the  bottom. 


Fig.  1  -  The  architecture  of  the  TEG  exhaust  heat  recovery 
system.  (1,4  -  Aluminum  block;  2  -  TEG  module  assembly; 
3  -  Compact  heat  exchanger  housing;  5  -  Diffuser). 


Coolant  water  circulates  through  the  flow  fields  inside  the 
aluminum  blocks  to  establish  the  temperature  difference.  TEG 
modules  are  assembled  between  the  housing  and  the  blocks. 

Apart  from  part  5,  the  other  four  parts  in  Fig.  1  form  the 
main  body  of  the  system.  The  size  of  the  main  body  is  deter¬ 
mined  by  the  number  of  TEG  modules  along  and  perpendic¬ 
ular  to  the  flow  direction,  namely  nrun  and  ncro.  These  four 
parts  are  discretized  into  segments  using  a  control  volume 
methodology.  As  illustrated  in  Fig.  2,  the  segments  are  shaped 
via  evenly  dividing  the  whole  system  along  and  perpendicular 
to  the  flow  direction  by  designated  numbers,  nx  and  ny.  Each 
segment  includes  a  part  of  TEG  modules  and  the  heat  reser¬ 
voirs  on  both  sides. 

Part  5  in  Fig.  1  is  the  diffuser,  the  function  of  which  is  to 
distribute  the  exhaust  gas  homogeneously  among  the  flow 
passages  inside  the  heat  exchanger  housing.  The  exhaust  gas 
is  from  a  1  kW  HTPEMFC  stack  operating  under  normal 
conditions.  The  gas  flow  rate  is  mgas  =  12.61  g/s.  Its  tempera¬ 
ture  is  tgas  =  148.2  °C.  Its  humidity  ratio  is  wgas  =  0.02076. 
Considering  that  the  stoichiometry  is  as  high  as  20  and  the 
oxygen  depletion  is  small,  it  is  of  satisfactory  accuracy  to 
model  the  exhaust  gas  as  ideal  moist  air.  The  thermodynamic 
data  are  from  Hyland  and  Wexler  [30].  The  remaining 
assumptions  of  the  numerical  model  are  as  follows: 

1)  There  is  no  heat  loss  from  the  system  to  the  surroundings. 

2)  The  electrical/thermal  contact  resistances  between  the  TEG 
modules  and  the  electric  load/heat  reservoirs  are  neglected. 
There  are  no  heat  bypasses  through  the  gaps  between  the 
TEG  modules. 
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3)  Thermoelectric  materials  distribute  homogeneously  inside 
each  module. 

4)  No  phase  change  occurs  inside  the  heat  exchanger.  Axial 
heat  conduction  is  neglected.  System  performance  is  iden¬ 
tical  in  the  direction  perpendicular  to  the  exhaust  gas  flux. 

5)  When  calculating  the  pressure  drop,  contraction  and 
expansion  pressure  losses  are  neglected. 

6)  Heat  conduction  and  convection  in  the  aluminum  blocks  as 
well  as  between  the  TEG  modules  and  the  aluminum  blocks 
are  not  taken  into  account.  The  coolant  water  flow  is  set  to 
be  high  enough  so  a  uniform  temperature  of  20  °C  can 
represent  the  cold  side. 

Based  on  the  above  assumptions,  the  governing  equations 
for  each  segment  can  be  given  based  on  the  energy  balance. 

2.1.  TEG  modules 

The  way  to  model  TEG  modules  is  reconsidered  compared  to 
other  publications  in  the  field  to  better  serve  the  system  level 
study.  Modules  are  treated  as  black  boxes  in  the  model.  Inside 
each  module,  thermoelectric  materials  are  supposed  homo¬ 
geneously  distributed.  Unlike  most  of  the  above  literature, 
module  parameters  are  from  the  experiments  in  [26],  not 
directly  substituted  by  material  properties.  Module  perfor¬ 
mance  calculated  from  empirical  parameters  is  more  repro¬ 
ducible  in  real-life  applications.  After  the  module  is 
discretized  in  the  x  and  y  directions  but  not  along  the  z  axis, 
TEG  performance  is  assumed  uniform  in  each  segment  but 
vary  with  temperature  between  segments.  For  system  design 
and  fabrication,  the  advantage  of  modeling  TEG  modules  in 
this  way  would  be  twofold:  a)  mitigating  the  model  complexity 
and  b)  keeping  the  model  precision  from  being  affected  by 
simplifying  the  details  inside  the  TEG  modules.  The  module 
parameters  used  in  the  current  system  are  listed  in  Table  1. 
The  Thomson  effect  is  neglected  because  the  temperature 
difference  is  relatively  moderate  in  this  case  [31]. 

For  each  TEG  segment,  the  corresponding  energy  balance  is 
described  by  the  following  equations: 


^h,TE  i 

J)  =  (Wi)  -  Tc,TE(i))/Rt,i  +  «iTh>TE(i)l,  -  0.5I?R€  j, 

(1) 

^c,TE  ( 

i)  =  (Th,TE(0  -  Tc,TE(i))/R,,i  +  aiTc,TE(i)li  +  0.5I?Re,i, 

(2) 

MO 

=  MtE  (0  _  Jc.TE  (0  • 

(3) 

All  segments  can  either  electrically  work  in  series  or  in 
standalone  mode.  In  standalone  mode,  the  external  ohmic 


resistance  equals  the  internal  when  one  segment  reaches  its 
peak  power  output  and  the  corresponding  electric  current  is: 

h  =  0.5 («i)(Th)TE(i)  -  TC)TE (i))/Re,i-  (4) 

The  total  power  output  of  the  whole  TEG  assembly  can  be 
summed  up  as: 

nx  ny 

Ptea  =  ]L£u;(0-  (5) 

2.2.  Compact  plate-fin  heat  exchanger 

Choosing  the  best  heat  exchanger  for  heat  recovery  is 
a  common  problem.  In  the  current  model,  a  lookup  table  is 
integrated  to  alleviate  this  problem.  The  table  contains  59 
types  of  commercialized  compact  plate-fin  heat  exchangers. 
All  the  heat  exchanger  geometry  and  relevant  parameters  in 
the  table  are  taken  from  [32].  Attempts  towards  numerical 
prediction  of  heat  exchanger  performance  easily  fail  in 
matching  experimental  data.  However  the  correlations  and 
parameter  data  herein  were  generated  from  experimental 
work  by  Kays  and  London.  These  correlations  and  data  have 
found  extensive  application  in  industry,  particularly  in  less- 
critical  designs.  Compared  to  homemade  heat  exchanging 
structures,  modeling  in  this  way  can  better  fulfill  the  system 
optimization  and  fabrication  requests. 

The  governing  equations  to  describe  the  heat  transfer 
between  the  TEG  assembly  and  the  hot/cold  fluids  in  each 
segment,  respectively,  are: 


^h,TE  (i 

)  —  ^gas  (0  —  Mas  Mass  ^gas  i  l)  • 

(6) 

^gas  (0 

=  £ctf  (0 Mas  (hgas.i  —  ^g*s,i)  > 

(7) 

^h,TE  (] 

i)=UAhx(i)(Tgas(i)-Th,TE(i)), 

(8) 

Tgas  ({ 

)  =  (Tgas(0  +  Tgas(i  + 1))/2. 

The  coefficients  in  the  above  functions  are  given  as  follows: 

l/UAtot(i)  =  l/UAhx(i)  +  l/UATE(i),  (9) 

£ctf  (i)  =  1  -  exp(-NTU(i)), 

NTU(i)  =  2UAtot  (i)  Tgas  (i)  /  mgas  (hgaS);  +  kgas,i+i)- 

hgas  is  the  enthalpy  of  the  exhaust  gas  and  h™  is  the  exhaust 
gas  enthalpy  at  the  temperature  of  the  coolant  water.  They  are 
calculated  on  each  segment  from  functions  found  in  [30]. 
UAhx(i)  is  the  heat  exchanger  conductance  in  each  segment, 
which  describes  both  the  air-exchanger  heat  convection  and 
the  heat  conduction  in  the  exchanger  body.  It  is  determined  by 
the  exchanger  material  and  structure,  the  local  fluid  proper¬ 
ties  as  well  as  the  flow  condition  of  the  exhaust  gas.  In  details, 
the  key  functions  are  [32]: 


Ahx(i)  =  z(i)HhxWhxdx/ny 

(ii) 

Co(i)  =  jhx^gas  (kgasj  +  kgas,i+i) /2Prj  ^  ^Tgas(i) 

(12) 

£ 

TO 

II 

TO2’ 

1 

(13) 

REj  —  ugasDh/^gas,i 

(14) 
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For  any  kind  of  heat  exchanger,  there  is  an  operational 
region  of  Reynolds  number.  This  information  is  also  available 
in  the  lookup  table. 

In  addition,  the  pressure  drop  through  the  heat  exchanger 
housing  can  be  evaluated  from  the  following  function, 

Aphx  ^gas/ ^Pgas,in^  ^/hxkhxPgas,in/  ^f^hPgas^ 

O'  )  ^  ^Pgas,in/ P  gas, out  ^  ^  )]  (15) 

The  Colburn  factor  jhx  and  the  Friction  factor  fhx  of  plate-fin 
surfaces  are  also  from  the  experimental  source  [32]. 

2.3.  Solution  methodology 

The  solution  methodology  of  the  numerical  model  is  imple¬ 
mented  by  using  a  discretized  e-NTU  approach  [33].  Input 
parameters  include  the  thermo-physical  properties  and  elec¬ 
trical  connection  style  of  TEG  modules,  configuration  param¬ 
eters  of  TEG  modules  and  plate-fin  heat  exchangers  as  well  as 
the  inlet  exhaust  gas  parameters  into  the  plate-fin  heat 
exchanger.  Then  the  fluid  properties,  thermoelectric  perfor¬ 
mance  and  heat  transfer  process  are  updated  on  each 
segment  along  the  flow  direction.  Finally,  results  are  exported. 

The  discretized  e-NTU  harnessed  here  is  to  make  sure  the 
model  is  better  suited  for  both  simulation-type  and  design- 
type  problems  [34].  Algebraically,  it  is  identical  with  the  log- 
mean  temperature  difference  method  (LMTD)  [34]. 


3.  Simulation  results  and  discussion 

In  the  following  numerical  simulation,  the  heat  exchanger 
type  used  is  called  ‘Strip-fin  plate-fin,  surface  l/4(s)-ll.l’  [35]. 
It  is  characterized  by  a  low  pressure  drop  for  high  heat 
exchanging  capability  [27].  More  details  of  its  specifications 
can  be  found  in  [32].  The  only  exception  is  in  Section  3.5, 
where  the  performance  of  different  heat  exchangers  is 
compared. 

3.1.  Model  validation 

A  check  is  carried  out  to  balance  the  computational  require¬ 
ments  and  the  result  accuracy.  There  are  10  Melcor  HT8 
modules  along  and  perpendicular  to  the  flow  direction,  under 
which  case  the  system  is  oversized.  All  TEG  segments  during 
the  simulation  run  are  in  standalone  mode.  In  other  words, 
each  segment  is  assumed  to  be  connected  to  a  matched  load 
to  maximize  its  power  output.  According  to  the  4th  assump¬ 
tion  above,  only  the  number  of  segments  along  the  flow 
direction  can  potentially  affect  the  simulation  results. 
Therefore  the  segment  number  per  TEG  module  is  varied  from 
1  to  20  to  observe  its  influence.  The  representative  results  of 
the  TEG  hot  side  temperature  and  the  system  power  output 
are  shown  below. 

In  Fig.  3,  it  is  clear  that  the  effect  of  the  number  of  segments 
on  the  TEG  hot  side  temperature  distribution  is  not  significant. 
If  this  temperature  distribution  was  studied  independently, 
there  would  be  no  need  to  pay  attention  on  the  number  of 
segments.  On  the  other  hand,  it  impacts  the  system  power 


output  evidently.  Fig.  4  illustrates  the  system  power  output 
and  its  relative  error  versus  the  segment  number.  For  a  devi¬ 
ation  less  than  0.5%,  cutting  each  module  into  6  segments  in 
the  flow  direction  is  the  minimum  requirement.  It  means  it  is 
not  necessary  to  discretize  each  module  by  the  width  of  each 
thermoelement  or  the  distance  between  adjacent  thermoele¬ 
ments  as  in  literature.  In  the  following  simulations,  this 
number  is  adjusted  accordingly. 

Simulations  are  then  conducted  to  validate  the  system 
model  with  experimental  data.  Firstly,  the  discretized  TEG 
model  is  compared  with  experimental  results  from  [26]. 
Throughout  the  validation,  three  temperatures  are  respec¬ 
tively  set  to  the  module  hot  side,  while  the  cold  side  is  kept 
constant  at  50  °C,  to  follow  the  experiment  settings.  Sorted  by 
the  temperature  difference,  the  module  power  outputs  are 
plotted  in  Fig.  5.  The  results  show  a  perfect  correspondence 
between  the  discretized  TEG  model  and  the  experiments. 

Secondly,  the  performance  of  the  heat  exchanger  part  is 
verified  against  the  experimental  results  from  [27].  The  heat 
exchanger  type  harnessed  in  the  experiments  was  also  the 
one,  ‘Strip-fin  plate-fin,  surface  l/4(s)-ll.l\  For  comparison, 
the  heat  transferred  and  the  pressure  drop  are  used  for 
comparison.  Results  are  shown  in  Fig.  6. 

According  to  Fig.  6,  it  can  be  concluded  that  the  simulation 
results  are  in  good  agreement  with  the  experimental  results. 
They  also  show  similar  trends.  The  difference  relatively 
observable  only  lies  in  the  low  flow  rate  region.  The  reason  is 
that  some  vital  details  are  missing  in  the  reference,  such  as 
the  exact  temperature  difference,  the  hot  gas  composition, 
among  others.  As  a  whole,  the  deviation  still  falls  inside  a  10% 
range.  So  it  is  concluded  that  the  accuracy  of  the  heat 
exchanger  model  is  acceptable. 

The  numerical  model  of  the  TEG  exhaust  heat  recovery 
system  model  is  thus  considered  developed  and  validated. 
Before  conducting  optimization  on  the  system  configuration, 
a  sensitivity  analysis  is  essential  to  identify  which  parameters 
have  more  significant  impact  on  the  system  performance. 

3.2.  Sensitivity  analysis 


A  multiple  parameter  sensitivity  analysis  is  carried  out  in  this 
section  to  determine  the  main  variables  of  this  system  that 


Fig.  3  -  TEG  hot  side  temperature  distribution  vs.  the 
segment  number. 
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Fig.  4  -  System  power  output  with  segments  from  1  to  20 
per  module. 


Hot  flux  flow  rate  (kg/s) 

Fig.  6  -  Heat  exchanger  model  validation  by  heat 
transferred  and  pressure  drop. 


affect  the  performance.  An  uncertainty  propagation  tool  is 
used  to  study  the  effects  on  system  model  output  when  the 
independent  parameters  vary  in  given  ranges  [35].  This 
method  was  initially  designed  for  experiments  to  assess  the 
deviations  in  measurement  of  the  calculated  variables.  As 
described  in  [36],  this  calculation  assumes  the  independent 
variables  are  uncorrelated  and  random,  i.e.  neglecting 
covariance  between  them;  the  uncertainty  in  the  system 
performance  can  be  expressed  as 

u,  =  ^(dv/dXi)2^,  where  v  =f(X1,X2,  (16) 

Here  factors  that  affect  the  system  power  output  PTEA.max 
and  the  heat  exchanger  pressure  drop  Aphx  are  analyzed  using 
the  uncertainty  propagation.  These  factors  can  be  distin¬ 
guished  into  three  categories:  system  structural  parameters, 
TEG  properties  and  operating  conditions.  System  structural 
parameters  further  have  two  subcategories,  namely  system 
dimensions  and  heat  exchanger  type.  Results  are  listed  in 
Tables  2  and  3.  Every  input  variable  is  set  with  a  relative 
uncertainty  of  5.00%. 

The  results  show  that  the  system  power  output  and  the 
pressure  drop  are  affected  by  an  uncertainty  of  17.08%  and 
13.46%,  respectively.  The  influences  are  rather  significant, 
which  means  that  some  of  the  input  parameters  should  be 


0  2  4 

Current  (A) 

Fig.  5  -  Validation  of  the  discretized  TEG  model. 


carefully  treated  in  the  simulation  for  higher  model  accuracy 
as  well  as  in  practice  for  better  system  performance. 

By  observation  it  can  be  seen  that  the  larger  sources  of  the 
power  output  uncertainty  are  five  parameters.  These  are  the 
exhaust  gas  temperature,  mass  flow  rate,  specific  heat,  heat 
resistance  of  each  TEG  module  and  the  module’s  Seebeck 
coefficient.  As  verified  by  many  literatures,  the  exhaust  gas 
temperature  has  extremely  significant  influence  on  the  TEG 
power  output.  The  exhaust  gas  temperature,  together  with  the 
mass  flow  rate  and  the  specific  heat,  determines  the  capacity 
of  the  heat  source,  which  is  available  for  TEGs  to  convert  into 
electricity.  Similarly,  the  module  heat  resistance  and  the 
Seebeck  coefficient  dominate  the  conversion  efficiency  of  the 
whole  TEG  assembly.  If  any  of  the  above  parameters  becomes 
larger,  the  system  power  output  will  increase  and  vice  versa. 

On  the  other  hand  some  input  parameters  have  negligible 
effect  on  the  system  power  output.  These  include  the  exhaust 
gas  pressure  and  most  system  structural  parameters,  such  as 
the  number  of  modules  crossing  the  flow  direction,  number  of 
modules  along  the  flow  direction  and  so  on.  The  power  output 
insensitivity  to  these  structural  parameters  actually  indicates 
that  the  current  system  is  oversized.  This  insensitivity  is  also 
the  side  effect  of  all  TEG  segments  working  in  standalone 
mode  in  the  current  system.  The  structural  parameters  will  be 
more  thoroughly  analyzed  in  the  following  sections.  Another 
parameter  with  little  effect  is  the  exhaust  gas  pressure.  This 
can  have  a  positive  effect  on  the  system  electrical  efficiency 
because  little  pressure  boost  is  required  for  the  recovery 
system.  This  can  also  simplify  the  connection  between  the 
recovery  system  and  the  fuel  cell  stack,  because  the  feed  gas 
pressure  also  has  a  minor  effect  on  the  fuel  cell  performance. 
As  a  whole,  the  system  parasitic  losses  are  decreased. 

Parasitic  losses  and  system  pressure  drops  are  roughly 
proportional  in  the  current  system.  From  Table  3,  it  can  be 
concluded  that  almost  all  the  uncertainty  of  the  pressure  drop 
comes  out  of  the  system  structural  parameters  and  the  exhaust 
gas  magnitude.  Among  them,  the  exhaust  gas  flow  rate  and 
number  of  TEG  modules  perpendicular  to  the  flow  have  the 
most  significant  effect.  The  influence  of  the  exhaust  gas  pres¬ 
sure,  the  height  of  the  heat  exchanger  and  number  of  TEG 
modules  along  the  flow  direction  is  also  rather  significant. 
These  influences  have  been  clearly  described  by  the  above  Eq. 


INTERNATIONAL  JOURNAL  OF  HYDROGEN  ENERGY  37  (2  0  12)  8490-8498 


8495 


Table  2 

-  Sensitivity  study  on  system  power  output. 

Input  parameters  and  description 

Value 

Uncertainty  contribution  (%) 

Afin\A 

Fin  area/total  area 

0.756 

-0.00 

Heat  transfer  area/total  volume  (m2/m3) 

915.6 

1.35 

Hhx 

Height  of  heat  exchanger  (m) 

0.00635 

0.71 

finthk 

Thickness  of  fins  (not  applicable  to  pin-fin)  (m) 

0.000152 

0.00 

COjqxau 

The  exchanger  average  heat  conductivity  (W/m2-K) 

80.58 

1.30 

ftcro 

Number  of  modules  cross  the  flow  direction 

10 

0.01 

ft  run 

Number  of  modules  along  the  flow  direction 

10 

0.01 

Oi  TE 

Seebeck  coefficient  of  each  TEG  module  (V/K) 

0.05 

19.67 

RTE,e 

Electric  resistance  of  each  TEG  module  (Q) 

2 

-4.92 

F-TE,t 

Heat  resistance  of  each  TEG  module  (K/W) 

1.54 

10.71 

ft^gas 

The  exhaust  gas  mass  flow  rate  (kg/s) 

0.01261 

7.96 

Cp  gas 

Average  exhaust  gas  specific  heat  (J/kg-°C) 

1047 

7.96 

wgas 

Humidity  Ratio 

0.02076 

-0.00 

Pgas 

Gas  pressure  (Pa) 

1 

0.00 

tgas 

Gas  temperature  (°C) 

148.2 

44.51 

tCw 

Coolant  water  temperature  (°C) 

20 

-0.89 

Output 

Variable  description 

Value 

Uncertainty  (%) 

PTEA,max 

Maximum  power  output  of  the  system 

24.10  (W) 

±  4.12  (W)  (±17.08%) 

(15).  On  the  contrary,  the  effect  of  the  exhaust  gas  temperature 
and  its  specific  heat  is  relatively  moderate.  It  can  be  noticed 
that,  even  in  the  current  oversized  system,  the  system  struc¬ 
tural  parameters  still  have  dramatic  impact  on  the  system 
pressure  drop.  To  lower  the  pressure  drop,  as  well  as  maximize 
the  system  net  power  output,  the  system  dimensions  and  the 
heat  exchanger  types  need  a  more  detailed  analysis. 

3.3.  Influence  of  ncro 

The  number  of  TEG  modules  perpendicular  to  the  flow  direction 
is  varied  to  investigate  the  impact  of  the  heat  exchanger  width 
on  the  system  power  output  and  the  pressure  drop.  All  TEG 
segments  in  this  study  worked  in  standalone  mode.  The 
number  of  TEG  modules  along  the  flow  direction  was  desig¬ 
nated  to  30.  Under  this  condition  the  system  was  always  over¬ 
sized  along  the  flow  direction  for  every  ncro.  In  the  simulation, 


ncro  increased  from  the  minimum  value  of  2.  Under  the  above 
settings,  the  system  power  output  will  reach  its  upper  limit  in  all 
cases. 

The  simulation  results  are  shown  in  Fig.  7.  The  width  tends 
to  decrease  the  electric  power  generated  and  the  pressure 
drop  simultaneously.  The  pressure  drop  falls  sharply  down  in 
the  beginning.  But  after  ncro  reaches  8,  the  effect  is  rather 
moderate.  On  the  other  hand  the  system  power  output  drops 
gradually  all  the  way  down  with  ncro  increasing,  although 
a  little  bit  faster  in  the  first  half  range.  The  increasing  width 
increases  the  cross  section  area  of  the  flow  passage.  The  flow 
speed  then  drops.  In  turn,  it  can  be  calculated  by  Eq.  (15)  that 
the  pressure  drop  will  decrease.  But  with  decreasing  flow 
speed,  Reynolds  number  is  also  becoming  smaller.  As  a  result, 
poorer  convective  heat  transfer  in  the  heat  exchanger  will 
reduce  the  temperature  difference  on  TEG  modules  and 
reduce  the  power  output. 


Table  3 

-  Sensitivity  study  on  the  system  pressure  drop. 

Input  parameters  and  description 

Value 

Uncertainty  contribution  (%) 

Afin\A 

Fin  area/total  area 

0.756 

0.00 

Heat  transfer  area/total  volume  (m2/m3) 

915.6 

-0.00 

Hhx 

Height  of  heat  exchanger  (m) 

0.00635 

-15.64 

finthk 

Thickness  of  fins  (not  applicable  to  pin-fin)  (m) 

0.000152 

-0.00 

COhXjau 

The  exchanger  average  heat  conductivity  (W/m2-K) 

80.58 

-0.00 

ftcro 

Number  of  modules  cross  the  flow  direction 

10 

-28.11 

ft  run 

Number  of  modules  along  the  flow  direction 

10 

13.13 

a!TE 

Seebeck  coefficient  of  each  TEG  module  (V/K) 

0.05 

0.00 

RrE,e 

Electric  resistance  of  each  TEG  module  (Q) 

2 

0.00 

F-TE,t 

Heat  resistance  of  each  TEG  module  (K/W) 

1.54 

0.01 

ft^gas 

The  exhaust  gas  mass  flow  rate  (kg/s) 

0.01261 

28.06 

Cp  gas 

Average  exhaust  gas  specific  heat  (J/kg-°C) 

1047 

0.03 

wgas 

Humidity  Ratio 

0.02076 

-0.00 

Pgas 

Gas  pressure  (Pa) 

1 

-13.80 

tgas 

Gas  temperature  (°C) 

148.2 

1.20 

tew 

Coolant  water  temperature  (°C) 

20 

0.02 

Output 

Variable  description 

Value 

Uncertainty  (%) 

A  Phx 

Pressure  drop  of  the  system 

383.10  (Pa) 

±  51.58  (Pa)  (±13.46%) 

8496 


INTERNATIONAL  JOURNAL  OF  HYDROGEN  ENERGY  37  (2  0  12)  8490-8498 


Number  of  modules  perpendicular  to  the  flow 


Fig.  7  —  Power  output  and  pressure  drop  vs.  number  of  TEG 
modules  cross  the  gas  flow. 


As  mentioned  above,  the  pressure  drop  accounts  for  most 
of  the  system  parasitic  losses,  through  determining  the  power 
demand  of  the  blower  in  the  system.  If  only  the  following 
function  from  [37]  is  used  to  assess  the  power  demand,  an  ncro 
value  higher  than  8  will  be  the  preferred  choice.  But  actually, 
the  blower  in  the  whole  fuel  cell  and  TEG  integrated  system 
mainly  serves  the  fuel  cell  stack.  Therefore  a  further  system¬ 
atic  study  is  essential  here  to  more  properly  predict  the  effect 
of  the  recovery  system  pressure  drop  on  the  net  power  output 
of  the  whole  system.  In  this  paper,  the  value  of  ncro  is  set  to  6  in 
the  following  analysis. 

Pcomp  —  jhgasBgasTinln(pout/Pin) //Mgas77comp  (17) 

3.4.  Influence  of  nrun 

A  simulation  was  also  carried  out  to  study  the  effect  on  the 
system  performance  from  the  length  of  the  flow  channels, 
which  is  determined  by  the  number  of  TEG  modules  along  the 
exhaust  gas  flow  direction,  nrun.  The  other  settings  were:  ncro 
was  fixed  to  6;  all  TEG  segments  were  either  in  series  or  in 
standalone  mode.  The  results  are  compared  in  Fig.  8. 

It  can  be  seen  that  the  two  system  power  output  curves  are 
overlapped  and  increase  very  fast  in  the  beginning.  Then  the 
power  output  curve,  when  all  TEG  segments  are  operating  in 
standalone  mode,  tends  to  approach  the  theoretical  upper 
limit.  There  is  no  optimum  value  of  nrun.  But  the  power  output 
of  all  segments  in  series  starts  to  drop  down  slowly  after 
reaching  the  maximum  point,  where  nrun  =  7.  This  effect  is 
similar  as  what  has  been  reported  in  [38].  It  was  believed  that 
the  peak  point  is  the  trade-off  between  the  electromotive  force 
and  the  internal  resistance,  when  all  TEG  modules  are  in 
series.  In  real-life  applications,  both  of  the  above  connection 
methods  of  the  TEG  segments  are  either  lacking  feasibility  or 
reliability.  A  mixed  connection  method  and  an  intermediate 
power  output  are  considered  more  feasible. 

On  the  other  hand  the  system  pressure  drop  under  both 
connection  conditions  is  roughly  proportional  to  nrun  in  the 
whole  range.  It  can  also  be  noticed  that  the  two  curves  are 
completely  overlapping.  This  phenomenon  corresponds  with 
the  above  sensitivity  analysis,  which  proved  that  the  system 
pressure  drop  is  independent  from  the  performance  of  the 
TEG  modules. 


Fig.  8  -  Power  output  and  pressure  drop  vs.  nrun. 


3.5.  Influence  of  heat  exchanger  design 

In  this  section,  different  types  of  heat  exchangers  are 
compared  by  the  influence  on  the  system  power  output  and  the 
pressure  drop.  Based  on  the  analyses  above,  6  TEG  modules 
were  aligned  across  the  exhaust  gas  flow  in  the  simulation;  7 
modules  were  along  the  flow  direction.  All  TEG  segments  were 
electrically  connected  in  standalone  mode.  Judged  by  their 
feasible  Reynolds  number  interval,  the  results  from  33  types  of 
plate-fin  heat  exchangers  among  all  59  types  are  shown  in 
Fig.  9.  The  results  are  sorted  by  the  power  output.  Heat 
exchangers  with  pressure  drop  larger  than  900  Pa  are  not 
included. 

It  can  be  seen  that  with  the  ascending  power  output,  the 
pressure  drop  fluctuates  significantly  between  different  heat 
exchangers.  Obviously,  heat  exchangers  with  higher  power 
output  and  lower  pressure  drop  are  superior.  From  Fig.  9,  four 
types  are  identified  of  interest  for  further  study,  namely  ‘Plain 
plate-fin,  surface  15.08’  (17),  ‘Pin-fin  plate-fin,  surface  PF-10(F)’ 
(25),  ‘Strip-fin  plate-fin,  surface  l/6-12.18(D)’  (28)  as  well  as 
‘Pin-fin  plate-fin,  surface  PF-4(F)’  (33).  The  type  ‘Strip-fin  plate- 
fin,  surface  l/4(s)-ll.l’  currently  used  in  the  exhaust  heat 
recovery  system  is  the  type  (19)  in  Fig.  9. 


Fig.  9  -  Power  output  and  pressure  drop  using  different 
heat  exchangers. 
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4.  Conclusions 

A  numerical  model  of  a  TEG  heat  recovery  system  for  the 
exhaust  gas  from  a  HTPEMFC  stack  was  developed.  It  is  based 
on  a  finite-element  approach  with  more  precisely  described 
gas  properties  and  heat  transfer.  For  both  design-type  and 
simulation-type  problems,  the  versatility  of  the  model  is 
enhanced  by  the  integrated  database  of  59  types  of  plate-fin 
heat  exchangers  and  the  discretized  e-NTU  solution  method. 
In  summary,  the  model  can  help  optimize  the  system  design 
and  operation  with  higher  accuracy  and  flexibility.  In  turn, 
this  can  benefit  the  system  fabrication. 

A  sensitivity  analysis  was  then  carried  out  on  the  system 
model.  It  is  found  that  exhaust  gas  parameters  and  heat 
exchanger  structure  have  a  significant  effect  on  the  system 
power  output  and  the  pressure  drop.  Gas  properties  should  be 
precisely  described.  The  configuration  of  the  system  should  be 
carefully  optimized  to  maximize  its  performance. 

Based  on  the  model,  the  optimized  system  configuration  of 
the  heat  recovery  system  is  finally  depicted.  There  are  6  TEG 
modules  crossing  the  exhaust  gas  flow  and  7  modules  along  the 
flow.  In  the  system,  the  TEG  module  type  is  assumed  to  be  the 
Melcor  HT8  type.  They  are  all  electrically  connected  in  series.  The 
heat  exchanger,  type  ‘Strip-fin  plate-fin,  surface  l/4(s)-ll.T,  is 
implemented,  which  was  also  chosen  by  other  researchers  [27]. 

Screening  the  whole  database  of  different  heat  exchangers, 
it  turns  out  that  the  heat  exchanger  structure  currently  used 
in  the  system  is  probably  not  the  best  choice;  Four  other  types 
are  identified  with  superior  performance  and  worthy  of 
further  analysis. 

Future  work  is  needed;  Experiments  are  still  needed  to  verify 
the  model  accuracy.  The  heat  transfer  phenomena  in  the 
coolant  water  side  and  water  vapor  phase  change  in  the  exhaust 
gas  could  be  included  to  improve  the  model  robustness.  In  the 
sensitivity  analysis,  uncertainty  intervals  of  the  input  parame¬ 
ters  need  further  consideration  to  better  reflect  realistic  condi¬ 
tions.  In  comparing  different  heat  exchangers,  the  system 
dimensions  will  be  optimized  for  them  respectively.  More  types 
of  heat  exchangers  can  also  be  included  in  the  database. 
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Nomenclature 

Aff.  minimum  free  flow  area,  m2 

Ahx:  total  heat  transfer  area,  m2 


c:  gas  specific  heat,  J/kg-K 

Dh:  hydraulic  diameter,  m 

fhx:  Friction  factor 

h:  exhaust  gas  enthalpy,  kj/kg 

H:  component  height,  m 

I:  electric  current,  A 

jhx:  Colburn  factor 

L:  component  length,  m 

nx:  total  number  of  segments  along  the  flow 

ny:  total  number  of  segments  crossing  the  flow 

P:  power  output  or  consumption,  W 

Ap:  pressure  drop,  Pa 

q:  heat  transferred,  W 

T,  t:  temperature,  K,  °C 

ugas:  core  mass  velocity,  Kg/s-m2 

u;:  power  output  of  each  segment,  W 

W:  component  width,  m 

Greek  symbols 

Ect f  the  effectiveness 
v:  kinematic  viscosity,  m2/s 
p :  exhaust  gas  density,  kg/m3 
a:  minimum  free  flow  area/frontal  area 
t:  heat  transfer  area/exchanger  volume,  1/m 

Subscripts/Superscripts 

c,TE:  cold  side  of  the  TEG  module 
cw:  coolant  water 
e:  electrical  resistance 
gas:  exhaust  gas 

h,TE:  hot  side  of  the  TEG  module 
hx:  heat  exchanger 
i:  segment  index 
max:  maximum  value 
t:  thermal  resistance 
TEA:  whole  TEG  assembly 

Abbreviation 


Co:  exchanger  heat  conductivity,  W/m2-K 

N TU:  number  of  transfer  units 

Pr:  Prandtl  number 

RE:  Reynolds  number 

UA:  heat  exchanger  conductance,  W/K 

ZT:  nondimensional  figure-of-merit  of  module 


